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We discovered kHz QPOs in 80 archived RXTE observations from the peculiar low-mass 
X-ray binary (LMXB) Circinus X-l. In 11 cases these appear in pairs in the frequency 
qq i range of ~230 Hz to ~500 Hz for the upper kHz QPO and ~56 Hz to ~225 Hz for the 

' lower kHz QPO. Their correlation with each other, which is similar to that of frequencies 

of kHz QPO pairs in other LMXBs containing a neutron star, and their variation by 
, a factor two confirm that the central object is a neutron star. These are the lowest 

■ frequencies of kHz QPO pairs discovered so far and extend the above correlation over a 

f*^ ' frequency range of factor four. In this new frequency range the frequency difference of 

\J^) ' the two kHz QPOs increases monotonically by more than ~170 Hz with increasing kHz 

, QPO frequency, challenging theoretical models. 

o , 

1. Introduction 

o ' 



Quasi periodic-oscillations (QPOs) are thought to reflect the motion of matter in 
the curved spacetime around neutron stars and black holes. The high frequency 
QPOs have proven to be a good indicator of the nature of the compact object. 
If the LMXB includes a neutron star, pairs of kilohertz (kHz) QPOs often appear 
' simultaneously with frequencies of several hundred Hz. Their frequencies vary as the 

system evolves but their difference appears to reflect the spin frequency of the star 
or half of it, in case where this is known from X-ray pulsations. Their coherences 
and amplitudes range up to about 200 and 20% respectively. On the other hand, 
black holes only rarely show high frequency QPOsQ and then mostly one at a time; 
the frequency ratio of the higher to lower frequency QPO appears to be about 3:2 
and they have small amplitudes and coherences. See Ref. 8 for a review. 

The frequencies of the two ('lower' and 'upper') kHz QPOs follow a tight rela- 
tion for all neutron stars from which they have been observed. 5 Most models for the 
generation of the kHz QPO pairs relate at least one of them to an orbital motion in 
the accretion disc. In the sonic-point spin-resonance model, 3 the upper kHz QPO 
is generated by a bright spot on the stellar surface generated by gas from a clump 
orbiting at the sonic radius, whereas the lower kHz QPO is produced by illumina- 
tion of gas at the spin-resonance radius. In the relativistic precession model, 6 the 
upper and lower kHz QPO are identified respectively as the nodal and periastron 
precession of gas clumps at the inner edge of the disc. Relativistic resonance mod- 
els 2 suggest resonances between orbital frequencies in the disc, which would produce 



a The term kHz QPOs does not apply to them since they never reach 1 kHz. 
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Fig. 1. Left: the frequencies of the upper kHz QPO (top) and of the low-frequency QPO (bottom) 
against that of the lower kHz QPO for our results from Circinus X-l as also for other sources from 
Ref. 5. Right: The frequency separation of the kHz QPO pairs for the same sources as before, as 
well as the best fit for a straight line (with a slope of 0.6±0.1) and the relativistic precession model 
to the Circinus X-l data. The slight misplacement of the latter in Ref. 1 is corrected here. 

kHz QPO frequencies at a fixed ratio, similar to what is occasionally seen in black 
hole systems. 

Circinus X-l is an LMXB whose nature was long disputed. Although Type- 
I X-ray bursts were observed from the field of the object, 7 its X-ray, radio, and 
spectroscopic properties are more typical of a black hole. Here we summarize the 
discovery of kHz QPO pairs reported in Ref. 1. These show it is a neutron star. 

2. Observations 

We have analyzed ~2 Msec of archival data taken during 1996-2005 using the Rossi 
X-ray Timing Explorer in the 3-60 keV energy range. We divided the observations 
into segments of 16, 48, 128, or 256 sec, Fourier-transformed them, and then aver- 
aged each one to get power spectra. We made a precise estimate of the deadtime 
values of the detector by fitting the Zhang function 10 to observations where Circi- 
nus X-l does not seem to contribute, and used those to subtract the Poisson noise 
power from all observations. After renormalizing the power spectra to rms squared 
we fitted a multi-Lorentzian function for the characteristic frequencies. By plotting 
the frequencies of the fitted QPOs from each observation against each other as in 
Ref. 9, we found two groups of correlated frequencies extending over several hundred 
Hz. Based on the high frequencies and the over an order of magnitude frequency 
variability we identified those as lower and upper kHz QPOs. We found in total 80 
observations with kHz QPOs. In 11 cases both appeared simultaneously. All their 
characteristics can be found in Ref. 1 ; in general they have relatively low coherences 
and amplitudes. 

In the left panel of Fig. [1] we plot the two frequencies. They follow the same 
relation that was found for other neutron stars in Ref. 5, but extend this relation to 
frequencies lower by a factor 4. It is of particular interest to see how the frequency 
separation behaves in this new frequency range for kHz QPOs. We see in the right 
panel of Fig.Q]that it increases over a factor 2 with increasing frequency of the lower 
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kHz QPO. This is different from what is seen in other sources that have shown 
kHz QPO pairs at higher frequencies. Mathematically speaking, the relativistic 
precession model and to a lesser extent the relativistic resonance model, can explain 
the observed behavior for a neutron star mass of 2.2±0.3M© or a 1:3 resonance 
respectiveljQ. The former is produced for slightly inclined, infinitesimally eccentric 
geodesic orbits around non-rotating compact objects, a framework that has been 
challenged. 4 Oscillations in the disk require an energy comparable to the binding 
energy of the star to produce relative amplitudes similar to those observed. The 
sonic-point spin-resonance model can explain the generation of kHz QPOs, but is 
challenged by the strongly increasing frequency separation with kHz QPO frequency 
and the low frequencies observed in Cir X-l. 



Acknowledgments 

This work was accomplished under the Marie-Curie Training Program HPMT-CT- 
2001-00245. S.B. acknowledges support by grants NSF AST 00-98399 and NASA 
NAG 5-12030 and the funds of the Fortner Endowed Chair at the University of 
Illinois. 



References 

1. Boutloukos, S., van der Klis, M., Altamirano, D., Klein-Wolt, M., Wijnands, R., 
Jonker, P., & Fender, R., Astroph. J. 653, 143 5 (2006); |astro-ph/ 0608089 

2. Kluzniak W. & Abramowicz M. (2003); |astroph/0304345| 

3. Lamb, F. K., & Miller, M. C, ApJL, submitted (2003); [astro-ph/0308179l 

4. Markovic, P., fc Lamb, F.K., Mon. Not. Ro. Astr. S., submitted (2000); 
|astro~ph/0009116p 

5. Psaltis, P., Bello ni, T., & van der Klis, M., Astroph. J. 520, 262 (1999); 
|astro-ph/9902013| 

6. Stella, L. & Vietri, M., Phys.Rev.L, 82, 17 (1999); |astro-ph/9803278"1 

7. Tennant, A. F., Mon. Not. Ro. Astr. S. 226, 971 (1987) 

8. van der Klis, M., in Compact Stellar X-ray Sources, Cambridge University Press, 39 
(2006) ; |astro-ph/041055T1 

9. van Straaten, S., van der Klis, M., Wijnands, R., Astroph. J. 619, 455 (2005); 
|astro-ph/0410505| 

10. Zhang, W., Jahoda, K., Swank, J. H., Morgan, E. H., Giles, A. B., Astroph. J. 449, 
930 (1995) 

11. Zhang, W., Morgan, E. H., Jahoda, K., Swank, J. H., Strohmayer, T. E., Jernigan, 
G., Klein, R. I., Astroph. J. 469, 29 (1996) 



b Note that the line for the relativistic precession model is slightly different than in Ref.l because 
of an error there in plotting the corresponding function. 



